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We investigate the magnetic field dependence of the physical properties of f^ -configuration 
systems with a crystalline-electric field (CEF) singlet ground state, which gives rise to a non- 
Fermi liquid (NFL) fixed point due to the competition between the Kondo-Yosida singlet and 
CEF singlet states. On the basis of the numerical renormalization group method, we find that 
the magnetic field breaks this NFL fixed point via two mechanisms: one causing the polarization 
of f-electrons and the other giving the "channel" anisotropy. These two mechanisms induce 
a difference in the magnetic field dependence of the characteristic temperature Tp(i7), the 
crossover temperature from NFL to Fermi-liquid behavior. While the polarization of f-electrons 
gives Tp(if) oc [x ^ 2.0), the "channel" anisotropy gives the if-independent Tp(if). These 
two mechanisms cross over continuously at approximately the crossover magnetic field He, 
where an anomalous iif-dependence of Tp(if) appears. Such T^(H) well reproduces the NFL 
behaviors observed in Thi_ccUccRu2Si2. We also find that the if-dependence of the resistivity 
and the magnetic susceptibility are in good agreement with the experimental results of this 
material. These results suggest that the NFL behaviors observed in Thi_ccUccRu2Si2 can be 
understood if this material is located in the CEF singlet side near the critical phase boundary 
between the two singlet states. 

KEYWORDS: Thi_ccUccRu2Si2, non-fermi liquid, magnetic field effect, f^-impurity problem, crystalline- 
electric field effect, Kendo effect, numerical renormalization group 



■ 1. Introduction 

In recent decades, non- Fermi- liquid (NFL) behaviors 
observed in some heavy fermion compounds and high- 
Tc cuprates have generated interest in issues on the 
] quantum-critical point (QCP). Of these NFL behaviors, 
" those based on a single correlated impurity in systems 
with the f^-configuration are classified into two sub- 
classes in which the QCP is triggered by the local critical- 
. ity: one is caused by the two-channel Kondo (TCK) effect 

■ due to the non-Kramers doublet state, and the other 
is caused by the competition between the crystalHne- 
electric field (CEF) singlet and the Kondo-Yosida (K- 
Y) singlet states. ^^^^^ Each of these mechanisms shows 
NFL behaviors below its characteristic energy scale 
because the systems flow toward an unstable fixed point. 
However, in real systems, small but relevant perturba- 
tions, leading the systems away from an unstable fixed 
point, give rise to a finite characteristic temperature Tp, 
the crossover temperature from NFL behavior to Fermi- 
liquid behavior. Namely, these two NFL behaviors are 
observed in the temperature (T) region < T < 
when T^<^T^. 

The NFL behaviors due to these two mechanisms 
are, in general, difficult to distinguish experimentally, 
especially in the case of U-impurity compounds, be- 
cause there exists some ambiguity in determining the 
CEF level scheme of U ions. Thi_^U^Ru2Si2 (x<0.07) 
is one such complicated heavy fermion impurity sys- 



tem. The NFL behaviors of this material are well scaled 
by impurity concentrations, so that many theoretical 
and experimental works have been carried out on the 
basis of these two mechanisms treating the U ion as 
an impurity.^' ^^^^^ In fact, the NFL behaviors of 
Thi_ajUajRu2Si2, such as the — InT divergence of both 
the magnetic susceptibility ximp and the Sommerfeld co- 
efficient 7imp = Cimp/T^ Cimp being the specific heat due 
to the impurity, and the anomalous temperature depen- 
dence of the resistivity Pimp, are consistent with those 
predicted by theories on the basis of these two mecha- 
nisms. In the case of Ri_a^Ua^Ru2Si2 (R=Y and La), the 
degrees of NFL behaviors are less prominent. Namely, 
Fermi liquid behaviors recover in the low temperature re- 
gions where Thi_a^Ua;Ru2Si2 exhibits prominent NFL be- 
haviors.^^' These differences can be understood from 
the viewpoint that the distances from the QCP are dif- 
ferent from compound to compound. In other words, 
Thi_a^Ua;Ru2Si2 is assumed to be accidentally located 
near the QCP. 

With the application of a magnetic field, however, 
there exist some aspects inconsistent with the NFL be- 
haviors on the basis of the TCK effect even in the case 
of Thi_a;Ua;Ru2Si2. First, if the NFL behaviors origi- 
nated from the TCK effect, the magnetic field would 
induce the increase in 7imp due to the release of the 
residual entropy by lifting the degeneracy due to the 

(2) 

doublet Fg ^ ground state of J = 4 orbitals in tetrago- 
nal symmetry. However, the suppression of the — logT 
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divergence of 7imp is observed in Thi_a^UajRu2Si2 by 
applying a magnetic field. Next, it was reported 
that Tp(iJ) of Tlii_a;Ua:Ru2Si2 shows an anomalous 
magnetic field {H) dependence, i.e., linear in in 
contrast to the quadratic dependence expected in the 
TCK model.^'^' '''^^ Considering these inconsistencies, 
it is troublesome to argue that the NFL behaviors in 
Thi_a^UicRu2Si2 can be explained by the theory based 
on the TCK effect.^^) 

In this paper, we study the magnetic field dependence 
of the NFL behaviors due to the competition between the 
K-Y singlet and the CEF singlet states in tetragonal sym- 
metry, and discuss its applicability to the magnetic prop- 
erties of Thi_ajUa:Ru2Si2. Yotsuhashi et al. have already 
discussed this problem on the basis of the two-orbital An- 
derson model with the "antiferromagnetic" Hund's rule 
coupling, the same as in the present paper. They have 
shown that the logarithmic increase in 7imp due to the 
competition between the two singlet states is suppressed 
by applying the magnetic field in a wide set of param- 
eters near the unstable fixed point, which is consistent 
with the experimental results of Thi_a;Ua;Ru2Si2 in a 
wide-temperature region. Here, we also take the same 
CEF scheme as that in ref. 11, and investigate the in- 
dependence of the magnetic susceptibility Ximp, the re- 
sistivity Pimp, and the characteristic temperature T^{H) 
obtained from these physical quantities. On the basis 
of the Wilson numerical renormalization group (NRG) 
method, we show that the iJ-dependence of Tp(iJ) 
changes at around the crossover magnetic field He, and 
that Tp (H) dit H ^ He reproduces the anomalous behav- 
ior observed in Thi_a^Ua^Ru2Si2. Namely, the anomalous 
properties in Thi_^U^Ru2Si2 can be fully explained by 
the present model. Moreover, the anomalous properties 
in Ri_a;Ua;Ru2Si2 (R=Y and La) are also consistent with 
our results. 

This paper is organized as follows. In §2, we introduce 
the model Hamiltonian to discuss the competition be- 
tween the two singlet states. In §3, the numerical result 
obtained by the NRG calculation of the magnetic field 
effect on Tp(i7), the resistivity pimp, and the magnetic 
susceptibility Ximp, are given in the cases of both singlet 
ground states. In §4, we discuss a scaling property of the 
iJ-dependence of T^{H) and its origin on the basis of 
the similarity of the unstable fixed point to the case of 
the TCK effect. In §5, we discuss the applicability of this 
scenario to the experimental result of Thi_a^Ua;Ru2Si2, 
and summarize our results in §6. 

2. Model Hamiltonian 

To discuss the competition between the two singlet 
states, we rewrite f^-states in the j-j coupling scheme us- 
ing f^-states in the j = 5/2 manifold. Here, we restrict the 
Hilbert space of f^ -states to two low- lying Kramers pairs 
and allot them the pseudospin states as follows: 




Fig. 1. CEF level scheme of low- lying states and their eigen- 
states. 

|r6,+) = | + i)^|o,t), (3) 

\Te,-) = \-l)^\OA). (4) 

The f^-states are also restricted to the four low- lying CEF 
states in the J = 4 manifold of tetragonal symmetry, 
which are written in the j-j coupling scheme within a 
manifold of j = 5/2 in f^ -configuration as follows: 



|r4) 




(I;,t>-It4», (5) 


|r3> 




(lt4) + i;,t)), (6) 


|ri!l) 


= /3|+3)-a|-l) = |t,t>, 


(7) 




= /3|-3)-a| + l) = |4,;). 


(8) 



Here, we assume the r4 singlet ground state shown in 
Fig. 1, where K and A represent the excitation energies. 
With the use of the pseudospin states (l)-(4), the f^-level 
scheme is reproduced by the "antiferromagnetic" Hund's 
rule coupling^ ^ ' 



HHund = -Y [StS^ + S^S+] + J.SfSI, (9) 

where coupling constants are defined sls J± = K and 
Jz = 2A — K. Sm is a pseudospin operator of the f- 
electron in the Hilbert space of the f^-state spanned by 
the orbitals m = 1 (Fl^^^) or 2 (Fg), and is defined as 

where fmcr is an annihilation operator of the f-electron 
in orbital m}^'^^^ The use of the j — j coupling scheme 
for f^-states is not necessary, in principle, for solving 
the present problem. However, it makes the problem 
more tractable in calculations based on the Wilson NRG 
method. 

Thus, the Hamiltonian is given by the two-orbital 
impurity Anderson model with the "antiferromagnetic" 
Hund's rule coupling as:^^'^^^ 

'H = He + Hhyb + % + '^Hund, (11) 

E E^fc4„.^fema' (12) 

m=l,2 

m=l,2 jfc^ 



J. Phys. Soc. Jpn. 



Full Paper 



Author Name 3 



/ J -^frnfmafma H" / ^ ^rnf^^fm"\-fjji^fml') 
ma m 

(14) 

where c^^^ is the annihilation operator of a conduction 
electron with the wave vector k and the spin a hybridiz- 
ing with the f-electron in the orbital m with a strength 
Vm- Efm and Um are the energy level of the f-electron 
and an intra-orbital Coulomb repulsion in orbital m, re- 
spectively, and the other Coulomb repulsion terms, like 
inter-orbital interaction, are implicitly included in the 
"antiferromagnetic" Hund's rule coupling (9). 

We consider the case when the magnetic field is ap- 
plied in the 2:- direction, the c-axis of Thi_2^Ua;Ru2Si2. 
The effect of the magnetic field for f^ -states is taken into 
account through the Zeeman terms defined by 

eeman 

m 

where the ^-factors of orbitals 1 and 2 are gi = 
90/49 and ^2 = 6/7, respectively. The effects of 
the magnetic field for f^-states (n = 2,3,4) are cal- 
culated using 'Hzeeman(f'^), which is the sum of the 
Zeeman term (15) for each f-electron. For example, 

(ri'i|Hzeema„(f')|r^'i) = ±(51+52)^8^^/2 for the f- 

state \TfX), i.e., (7) and (8), and {UA |^Zeema„(f^)| U 

,t) = 92mH/2 for the f^-state | UA) = flfflA-f\0), 
where |0) is the vacuum state, and so on. In the same 
manner, the Van Vleck contribution arising from the 
off-diagonal term between r4 and Fs in the f^-singlet 
manifold is estimated as (F3|'Hzeeman(f^)|F4) = —{gi — 
g2)/J^BH/2. However, this value is much smaller than that 
estimated in the LS'-coupling scheme in the J = 4 man- 
ifold, (F4I - gjmJzHlTs) = -2gjiiBH^ with gj = 4/5, 
because the higher F^^^ doublet state in the j = 5/2 
manifold in the f^ -configuration has been discarded in 
constructing our model Hamiltonian shown by eqs. (11)- 
(14). In fact, if we construct the F4 and F3 singlet states 
in the j = 5/2 manifold as 




the off-diagonal term is estimated as 

(r3|^Zeeman(f')|r4) = 2gj^BH with g^ = 6/7, 

which almost coincides with the value estimated in the 
I/S'-coupling scheme in the J = 4 manifold. Thus, to 
take into account the Van Vleck contribution properly, 
we adopt the off- diagonal matrix element in the J = 4 
manifold other than the contribution of the f^ -based 
Zeeman term. 

We transform the conduction band part of the Hamil- 



tonian (12), with a logarithmic discretization parame- 
ter, A = 2.5, into the one- dimensional semi-infinite chain 
model and carry out the Wilson NRG method. For 
simplicity, we take conduction bands to be symmetric in 
the energy space (with an extent from —D to D) centered 
at the Fermi level. We keep the low-lying 4000 states in 
each iteration step. 

3. Characteristic Temperature T^{H) 

The Hamiltonian (11) has two stable fixed points. One 
is the K-Y singlet fixed point (KY SFP) where the spin 
degree of freedom of each f-electron is screened by the 
conduction electrons with the same symmetry as the f- 
electron, leading to the phase shift in the unitarity limit 
as = 7r/2 (m = 1,2). The other is the CEF sin- 
glet fixed point (CEF SFP) where two f-electrons form 
the singlet state due to the CEF effect, characterized by 
= (m = 1,2). Along the boundary of these two 
stable-fixed-point regions, there exists a locus of the un- 
stable fixed points across which the ground state is in- 
terchanged. Around this line, NFL behaviors appear at 
T * < T < = min(TK2,i^), where Tk2 is the lower 
Kondo temperature of two f-orbitals. 

In general, Efm and Um, the energy level and the 
Coulomb interaction of each f-orbital, respectively, are 
different. However, for simplicity, we take the same val- 
ues for each orbital, and the difference in characters of 
each orbital is introduced only through Vm- The Kondo 
temperature of orbital 2 is postulated to always be lower 
than that of orbital 1, i.e., Tki > Tk2, and the param- 
eters of the Hamiltonian (11) are fixed as Efi = Ef2 = 
-0.4, Ui = U2 = 1.5, Vi = 0.45 and V2 = 0.30 in the unit 
of D throughout this paper. In addition, the magnetic 
field H is measured in the unit of D/ ji^- In the case of 
if = A = 0, the Hamiltonian (11) reduces to two in- 
dependent impurity Anderson models, where the Kondo 
temperatures determined by the definition of Wilson, i.e., 
4TKXimp(T = 0) = 0.413, are Tki = 4.52 x 10"^ and 
Tk2 = 3.43 X 10-3. In this paper, we set the CEF level 
splittings as if > A, which reproduces the anisotropy of 
the magnetic susceptibility, Xz > as pointed out in 
ref. 11. Moreover, we fix A = 0.12 and control the degree 
of the competition by varying the CEF level splitting K. 
For the parameter set above, 0.0464 gives an un- 

stable fixed point, i.e., the ground state is the K-Y singlet 
for K kK"" and the CEF singlet for K > K\ 

Figure 2 shows the temperature dependence of the spe- 
cific heat in two cases: K = 0.0440 in the KY SFP region 
and K = 0.0488 in the CEF SFP region. The character- 
istic temperature Tp is defined as the lowest temperature 
at which the specific heat Cimp(^^) = dSimp{H)/dlnT, 
where S'imp is the entropy due to the impurity, has a 
peak corresponding to the release of log y/2 entropy that 
characterizes the unstable fixed point. For these pa- 
rameters, the characteristic temperatures are obtained 
as - 1.44 X 10"^ in the KY SFP region and as 
T* - 1.33 X 10-^ in the CEF SFP region. In the case of 
the KY SFP region, Tp shghtly decreases for a magnetic 
field H = 1.0 X 10"^, but increases for the other values 
of the magnetic field. On the other hand, in the case of 
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Fig. 2. (Color online) Temperature dependence of the specific heat Cimp for a series of magnetic fields H (0 < H < 3 x IQ-^) for (a) 

1.33 X 10-6 in the CEF SFP. 



K = 0.0440 with T* = 1.44 x 10"^ in the KY SFP, and (b) K -- 
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Fig. 3. (Color online) Characteristic temperature Tp vs K = 
[K — K*)/K* for a series of magnetic fields H. Data points 
in the KY SFP are represented by closed symbols, while those in 
the CEF SFP are represented by open symbols. Arrows indicate 
the positions of the dip of Tp(if). 

the CEF SFP region, Tp increases at all values of the 
magnetic field. 

Figure 3 shows the K = {K - K*)/K* dependence of 
T^{H). The characteristic temperature Tp is decreased 
by the competition. One can see in Fig. 3 that the it- 
dependence of Tp at ili = is given by Tp ex K'^ around 
the QCP, indicating that Tp gives a degree of deviation 
from the QCP. When the magnetic field is applied, Tp 
increases, and the energy spectrum no longer suddenly 
interchanges at = if* because the ground state is 
the mixed state between the CEF singlet and K-Y sin- 
glet states. In the KY SFP region, there is a dip (indi- 
cated by arrow in Fig. 3) at which Tp takes a minimum 
but remains non-zero. The energy spectrum obtained by 
the NRG calculation "gradually" crosses over between 
the types of the CEF and the K-Y singlet states around 
this dip. Namely, it is the point where the dominant sin- 
glet state of the two singlet states interchanges. As the 
magnetic field increases, this dip moves from K = if* 
to the low K region, which indicates that the magnetic 
field increases the weight of the CEF singlet state com- 
pared with that of the K-Y singlet state. This increase in 
the weight of the CEF singlet state originates from the 



off-diagonal term between the Fa and r4 f^-CEF singlet 
states because it stabilizes the energy level of the r4 CEF 
singlet ground state. Hereafter, we investigate the in- 
dependence of physical quantities in two cases being close 
to the QCP: K = 0.0440 with T^ = 1.44 x 10"^ in the 
KY SFP region and K = 0.0488 with T^ = 1.33 x 10"^ 
in the CEF SFP region. 

Figure 4 shows the frequency dependence of the to- 
tal scattering rate l/r{uj) at T = in two cases, i.e., 
K = 0.0440 and K = 0.0488, where l/r{uj) is the sum of 
contributions from each orbital, and the spin component 
l/rma{^) = 27r\Vm\'^Ama{oj), ^ma(^) being the single- 
particle spectral function. Data points in the KY SFP 
region are represented by closed symbols, while those 
in the CEF SFP region are represented by open sym- 
bols unless stated explicitly. The T-dependence of the 
resistivity can be inferred from l/r{uj) because uj and 
T are of the same order in the Fermi-liquid theory, 
e.g., in the case of the single orbital Anderson model, 
l/Ticu, T) ^ [l/r(0, 0)] [1 - + t,^T^)/2,TI + • • • ] .^D 
Although it is not shown as figures in the present pa- 
per, 1/t{uj) increases logarithmically in the region of 
UJ > Tki owing to the screening of the f-electron by con- 
duction electrons in orbital 1 in both cases. With de- 
creasing UJ toward Tp(<C Tk2 < Tki), shows a 
logarithmic increase in the KY SFP region, but it shows 
a logarithmic decrease in the CEF SFP region. Finally, 
the Fermi- liquid behavior is restored in both regions, i.e., 
l/r{uj) (X cj^ at < T^{H = 0). 

When the magnetic field is applied, the residual scat- 
tering rate I/tq = l/r(a;)|c^^o decreases in the KY SFP 
region, as seen in Fig. 4(a). There are two origins that in- 
duce such an increase in 1 /tq : one is the mixing between 
the K-Y and CEF singlet states in the ground state, and 
the other is the mixing between the r4 and Fs singlet 
states through the off-diagonal term. As a result, the 
magnetic field leads to the polarization of each f-electron. 
These magnetic moments make a singlet state, leading to 
the reduction in the phase shift. On the other hand, in 
the CEF SFP region, I/tq increases because the CEF- 
type ground state is polarized, and its magnetic moment 
scatters off conduction electrons, leading to an increase 
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Fig. 4. (Color online) Frequency dependence of the total scatteri 
for (a) K = 0.0440 with T* = 1.44 x 10"^ in the KY SFP, and 

in the phase shift and l/r{uj) at < Tp(i^ = 0). Because 
the weight of the CEF singlet state in the ground state 
markedly increases compared with that of the K-Y sin- 
glet state around the QCP, l/r{uj) in the KY SFP region 
shows the same T-dependence as in the CEF SFP region 
under a high magnetic field. After all, the T-dependence 
of the resistivity pimp due to the impurity scattering is 
essentially given by that of l/r{uj = T). 

The magnetic susceptibilities Ximp = dM/dH are 
shown in Fig. 5 in these two cases. The magnetiza- 
tion M consists of Mi [arising from the Zeeman term 
'Wzeeman(f^)] and M2 [arising from the Van Vleck term 
in the f^ configuration, eqs. (5) and (6)]. Mi is given as 
the thermal average of the magnetic moment m, which is 
calculated as m = ±(^1 +^2)mb/2 for f^-state |r5^^), i.e., 
eqs. (7) and (8), and m = /iB^2/2 for f^-state | ti^ t)^ and 
so on. On the other hand, M2 is given by the effect of the 
off-diagonal element of the magnetization between the f^- 
CEF singlet states Fs, i.e., eq. (6), and i.e., eq. (5). In 
both cases, Ximp(T) shows the logarithmic T-dependence 
at approximately T ~ Tki and Tp < T < min(TK2,i^)- 
The magnetic field reduces the coefficient of the — log T 
term at Tp < T < min(TK2, K) and the Van Vleck con- 
tribution. In particular, these reductions in the KY SFP 
region are smaller than those in the CEF SFP region. 
The origin of this phenomenon is the interchange of the 
weight of the two singlet states in the ground state. At 
higher magnetic fields, in both cases, a broad peak ap- 
pears at T ~ Tp(i^), where Tp is obtained from Cimp, as 
mentioned above. 

4. Scaling Behavior of Characteristic Tempera- 
ture T*(ff) 

From the data of the T 01 uj dependence of Cimp, Ximp, 
and 1/r under the magnetic field we obtain the in- 
dependence of T^{H)^ as shown in Fig. 6(a), where the 
Tp{H) of 1/r and Ximp are defined as the temperature 
at which the logarithmic T-dependence stops with de- 
creasing uj and T toward 0. We define the crossover mag- 
netic field He as the intersection point of linear fits (dot- 
ted lines) on the log-log plot for high and low magnetic 
field regions in the CEF SFP region {K = 0.0488 and 
K = 0.05172), as shown in Fig. 6(a). It is remarkable 



rate l/r((ju) for a series of magnetic fields if(0<i^<3xl0 ^) 
K = 0.0488 with T* = 1.33 x 10"^ in the CEF SFP. 

that the thus-determined H^s for three different physi- 
cal quantities, i.e., Cimp, Ximp, and 1/r, almost coincide 
with each other, as shown by arrows in Fig. 6 (a), giving 
a solid basis for defining the crossover magnetic field He. 
In the KY SFP region {K = 0.0440 and K = -0.05172), 
the HcS are similarly defined, although Tp obtained from 
Ximp exhibits a tiny dip near H = He. Linear fits are not 
shown for presentation clarity. The crossover magnetic 
field He so determined for Cimp, 1/t, and Ximp almost 
coincide again. These H^s in the two singlet-fixed-point 
regions almost coincide with each other because the abso- 
lute values of K for these two parameters are almost the 
same, while He depends crucially on = {K — K*)/ K* ^ 
the deviation from the QCP. The normalized character- 
istic temperature, T^{H)/T^{H = 0), is shown in Figs. 
6(b) and 6(c) as a function of the normalized magnetic 
field H/ He from the three physical quantities mentioned 
above. Figure 6(b) shows the results for the KY SFP with 
parameters K = 0.0440 giving T^{H = 0) = 1.44 x 10"^ 
and K = 0.0460 giving T*(i7 = 0) = 2.80 x 10-^ and 
Fig. 6(c) shows those for the CEF SFP with parame- 
ters K = 0.0488 giving T^{H = 0) = 1.33 x 10"^ and 
K = 0.0468 giving T^{H = 0) = 3.79 x 10"^ These four 
lines exhibit a good scaling property in both cases, which 
indicates the following two important facts. 

First, the Tp(i7)'s of the three quantities Cimp, 1/t, 
and Ximp exhibit qualitatively the same behaviors, while 
they are qualitatively different. Second, the scaling prop- 
erty holds in both the CEF SFP and KY SFP regions 
even if the degree of the deviation from the QCP is 
different. Note that the shapes of the normalized plots 
for the two stable fixed points are different from each 
other. In the region H < He{K), T^{H) is indepen- 
dent of H and T^{H) is robust particularly against the 
low magnetic field, so that the three physical quantities 
discussed above are not affected appreciably. We previ- 
ously found this robustness of Tp(iJ) against H^ as re- 
ported in ref. 13. On the other hand, all the Tp(i^)'s 
show the i7-dependent form as T^{H) ex H^ in the re- 
gion He{K) <^ H < min(TK2,^). The exponent x is 
estimated to be x 2.0 in both the CEF and KY SFP 
regions. These two regions continuously cross over at ap- 
proximately H He{K). Of course, it is possible that 
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Fig. 5. (Color online) Temperature dependence of the susceptibilit 
K = 0.0440 with T* = 1.44 x 10-^ in the KY SFP and for (b) K 

T^{H) can be fitted as T^{H) ex i7 in a very narrow 
region of the magnetic field near H = H^^ especially in 
the CEF SFP region, as shown in Fig. 6(c). However, 
such a scaling behavior should be regarded as that of a 
crossover, but not an asymptotic anomalous behavior. 

Such a i^-dependence of T^{H) can be understood by 
considering the similarity of the unstable fixed point to 
that of the TCK effect. In the case of the TCK model, 
there are two origins that break the unstable fixed point, 
the magnetic field that polarizes the local spin leading 
to the "unusual" Fermi-liquid fixed point characterized 
by the energy scale Tp ex i^^/Tx,^^ and the channel 
anisotropy of the exchange interaction that leads to the 
Fermi-liquid fixed point.^'^'^^ Indeed, the Hamiltonian 
(11) can be regarded as the TCK model below Tki be- 
cause the f-electron in orbital 2 interacts with two "con- 
duction" electron channels: one is the conduction elec- 
trons in orbital 2, and the other is a complex of con- 
duction electrons and the f-electron in orbital 1 that is 
screened by the conduction electrons with the same sym- 
metry as the f-electron in orbital 1.^^'^^) In the present 
model, these two types of "conduction" electron serve as 
channels. The change in the energy difference between 
the two singlet states, which is induced by the mag- 
netic field, as mentioned above, affects the coupling con- 
stants between the f-electron in orbital 2 and the "con- 
duction" electrons. Namely, these two coupling constants 
exhibit magnetic field dependences. Thus, the "channel" 
anisotropy for the two types of "conduction" electron is 
induced by the magnetic field, and the system goes to the 
Fermi-liquid fixed point, even though the ground state is 
a mixture of the two singlet states. In other words, in 
the present model, the magnetic field breaks the unsta- 
ble fixed point via two mechanisms, the polarization of 
f-electrons and "channel" anisotropy. 

In the region H <^ He, the system flows into the Fermi- 
liquid fixed point induced by the "channel" anisotropy. In 
the CEF SFP region, Tp shows little change against a low 
H. However, in the KY SFP region, Tp slightly decreases 
as H increases, which corresponds to the dip in Fig. 3. 
Namely, the weight of the CEF singlet state in the ground 
state increases compared with that of the K-Y singlet 
state, and the dominant singlet state of the two singlet 



y Ximp for ci series of magnetic fields if (0 < if < 3 x 10 ^) for (a) 
= 0.0488 with T* = 1.33 x 10"^ in the CEF SFP. 

states interchanges H = H^,. On the other hand, in the 
region He <^ the magnetic field induces the Fermi- 
liquid fixed point by the polarization of f-electrons, as 
in the case of the TCK effect, because T^{H) is charac- 
terized by i^^, as in the case of the TCK effect under 
the magnetic field. These two effects compete with each 
other at approximately H ~ Hc{K), giving the crossover 
between the two regions. 

The exponent of H in Tp(i^) asymptotically ap- 
proaches 2.0 in the high magnetic field region H ^ He in 
both SFP regions, as shown in Figs. 6(b) and 6(c). How- 
ever, the magnetic field necessary to reach the Tp oc H^ 
behavior in the CEF SFP region is higher than that in 
the KY SFP region. This difference stems from the exis- 
tence of the Fs excited CEF singlet state, which gives an 
additional magnetic field dependence for the r4 CEF sin- 
glet ground state through the off-diagonal term between 
these two CEF singlet states. This is verified by a NRG 
calculation that the exponent of H in the CEF SFP re- 
gion readily comes close to 2.0 at > i^c, as in the KY 
SFP region, if we discard the off-diagonal term between 
f^-CEF singlet states by the magnetic field, although no 
explicit result is shown here. 

5. Comparison with Experiment on 
Thi_^U^Ru2Si2 

In the CEF SFP region near the QCP, l/r{uj) shows 
a behavior consistent with the results of Pimp(^) in 
Thi_a;Ua;Ru2Si2. Namely, l/r{uj) exhibits a logcj-like 
decrease toward Tp and increases as H increases corre- 
sponding to the positive magnetic resistance. Note that 
Tki is considered to be much higher than Tk2 because 
the log T-like increase in the resistivity at approximately 
T ^ Tki is not observed in this material. This result 
is consistent with that obtained by Yotsuhashi et al. who 
showed that l/r(ct;), Ximp(^), and the i^-dependence 
of 7imp reproduces these physical quantities observed in 
Thi_a^Ua;Ru2Si2. A more remarkable finding is that the 
JY-dependence of T^{H) at approximately H ~ Hc{K) 
for K = 0.0488 in the CEF SFP region reproduces the 
i^-dependence of T^{H) observed in Thi_a^Ua^Ru2Si2 an- 
alyzed from the resistivity. As shown in Figs. 7 and 
8, our theoretical results reproduce almost perfectly the 
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Fig. 6. (Color online) (a) Characteristic temperature Tp(if) vs 
magnetic field H in the two cases with CEF (for K = 0.0488 and 
K = 0.05172 shown by open symbols) and K-Y (for K = 0.0440 
and K = —0.05172 shown by closed symbols) singlet fixed points. 
The crossover magnetic fields HcS are shown by down arrows in 
the case of the CEF singlet fixed point, (b) and (c) are scaling 
plots of T*(if)/T*(/f = 0) vs H/Hc for the KY SEP and CEF 
SEP, respectively. 

experimental observation on the normalized magnetic 
field H/Hi dependence of T^{H) normalized by Tp(i^i), 
where Hi = 3.0 x 10~^D//j.b {D being half the band- 
width of conduction electrons) for our theoretical result 
and Hi = 1[T] for experimental result in ref. 18. Namely, 
our scahng plot is in good agreement with that of the ex- 
perimental result of Thi_a^Ua^Ru2Si2.^^^ 

It is also emphasized that our theoretical analysis 
strongly suggests that Thi_a;U^ERu2Si2 is located in the 
CEF singlet side near the critical phase boundary be- 
tween the KY SFP and CEF SEP regions. On the other 




Fig. 7. (Color online) Comparison between theoretical and ex- 
perimental results of T-^(H). The red o symbols are for the in- 
dependence of Tp(if) (upper and right scales) obtained theo- 
retically from Cimp for K = 0.0488. The blue □ symbols are 
for the if -dependence Tp (H) (lower and left scales) observed in 
Thi_xUa:Ru2Si2 for the resistivity pimp, which was scaled lin- 
early by Toth et al. in ref. 18. 
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Fig. 8. (Color onhne) ii/iii vs T* (ii)/T* (iii), normalized char- 
acteristic temperature, for theoretical and experimental results. 
Those obtained from Qmp and pimp are normalized by the value 
at iii = 3.0 X 10~^, while the experimental result is normalized 
by the value at iii ^ 1[T]. 

hand, it is emphasized in ref. 18 that Tp{H) is propor- 
tional to the magnetic field H^ especially in the analysis 
of the T-dependence of Ximp- However, it is apparent 
that the Hnear fit used in ref. 18 fails to reproduce the 
experimental results for the T^{H) obtained from the 
resistivity pimp- The statement in ref. 18 stemmed from 
the analysis of Ximp assuming that the coefficient of the 
— log T term in Ximp were independent of H. Experi- 
mentally, however, magnetic fields up to 5 Tesla seem 
to change this coefficient and markedly reduce the Van 
Vleck contribution to Ximp-^^^ 

In the case of Ri_ajUajRu2Si2 (R= La and Y), the 
experimental results show that the Tp's of these ma- 
terial are higher than that in the case of R=Th.^^'^^^ 
This indicates that parameter sets of these materials 
may be located more distant from the QCP than that 
of Thi_a^U^ERu2Si2. Our theoretical result predicts that 
pressure may induce the transition from the CEF SFP 
region to the KY SFP region, giving rise to a marked 
increase in I/tq. 

In the present paper, we take the same CEF level 
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scheme as that discussed in ref. 11, because such a 
level scheme can reproduce the experimental results in 
Ri-a;Ua;Ru2Si2 (R=Th, Y and La). However, even if the 

(2) 

low-lying CEF scheme is r4-r^"^-r3, we have verified 
that similar NFL behaviors and the magnetic field de- 
pendence of Tp(iJ) occur although the results are not 
shown in the present paper. This indicates that similar 
NFL behaviors would be obtained if there exist a CEF 
singlet ground state and a strong hybridization between 
conduction electrons and the f-electron, namely, details 
of the CEF scheme would not be essential matters. Note 
that CEF states with the T2 singlet ground state, pro- 
posed as a plausible candidate for the "Hidden Order" 
state of URu2Si2,^^'^^^ would exhibit the local non- Fermi 
liquid behaviors discussed in the present paper. The ac- 
tual calculation in those CEF level schemes is left for 
future study. 

Quite recently, it is argued that the iJ-dependence of 
NFL behaviors observed in Thi_a^Ua^Ru2Si2 can be un- 
derstood within the TCK model if some conditions would 
be satisfied in the CEF level scheme adopted in refs. 22 
and 23. A crucial difference between the result in ref. 
24 and our present one is the magnetic-field dependence 
of Tp in the high magnetic field region H > 5 Tesla. 
Namely, there exists a region of a magnetic field where 
Tp oc in the high magnetic field region in our result 
that is consistent with experiments as shown in Fig. 7, 
while there exists no such a region of a magnetic field in 
the scenario of ref. 24. 

6. Conclusions 

We have investigated the magnetic field {H) depen- 
dence of the NFL behaviors arising from the competi- 
tion between the CEF and K-Y singlet states in tetrag- 
onal symmetry. We have found that the characteristic 
temperature Tp(iJ), below which the Fermi liquid be- 
havior recovers, changes its iJ-dependence at approxi- 
mately the critical magnetic field He. While Tp is not 
affected by H in the region H ^ it is expressed as 
T^ (X H"" {x 0^ 2.0) in the region H > H^. For both 
high and low magnetic field regions, such a behavior of 
Tp(iy) follows the scaling form even if the degree of the 
deviation from the QCP is different. We have also found 
that Tp(i^) shows an anomalous iJ-dependence at ap- 
proximately H ^ He, which is in good agreement with 
that observed in Thi_^U;ERu2Si2. 
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